Revista Internacional de 

Ingenieria 
de Estructuras 


Vol.22, 2 135-162 (2017) 


MODELADO DE DISIPADORES HISTERETICOS DE 
ENERGIA EN ESTRUCTURAS DE CONCRETO 
REFORZADO. DIFERENCIAS ENTRE ANALISIS 
EN EL PLANO Y EN TRES DIMENSIONES 


Horacio de Jesus Nangullasmu Hernandez 

Candidato Doctoral, Posgrado en Ingenieria Estructural 

Universidad Autonoma Metropolitana-Azcapotzalco 

Av. San Pablo 180, Col. Reynosa Tamaulipas, 02200, Ciudad de Mexico 

e-mail: hnanqu@hotmail.com 

Arturo Tena Colunga 

Profesor, Departamento de Materiales, Universidad Autonoma Metropolitana-Azcapotzalco 
Edificio P4, 3er Piso, Av. San Pablo # 180, 02200 Ciudad de Mexico 
Correo electronico: atc@correo.azc.uam.mx . 


RESUMEN 

Se evaluan los resultados obtenidos entre dos programas de analisis estructural, a 
partir de analisis estaticos no lineales y dinamicos paso a paso, con el objetivo de 
evaluar el efecto de ortogonalidad en el modelado de marcos de concreto 
reforzado con disipadores histereticos de energia apoyados en contraventeo 
chevron. Se estudiaron marcos pianos y modelos en tres dimensiones con el 
principal objetivo de comparar sus resultados. Estos se valoran a partir de curvas 
de capacidad, mapeos de fluencias inelasticas, periodos de vibrar y ciclos de 
histeresis de los disipadores utilizados. Finalmente, se presentan comentarios 
acerca del modelado riguroso en tres dimensiones y del impacto de la rigidez del 
nudo viga-columna en la rigidez global de la estructura. 

Palabras clave: disipadores de energia, fusibles estructurales, estructuras de 
concreto reforzado, analisis estaticos no lineales, analisis dinamicos paso a paso 


MODELLING OF HYSTERETIC ENERGY DISSIPATION 
DEVICES IN REINFORCED CONCRETE STRUCTURES. 
DIFFERENCES BETWEEN 2D AND 3D ANALYSES 


ABSTRACT 

Results obtained from pushover and nonlinear time-history analyses are compared 
in order to evaluate orthogonality effects in reinforced concrete frames with 
hysteretic energy dissipation devices. Plane and three-dimensional frames were 
modeled using two different software. Their results were evaluated from capacity 
pushover curves, inelastic demands mappings, fundamental periods and hysteresis 
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loops of the dissipation devices. Finally, comments are presented about rigorous 
modeling in three dimensions related to the importance of the stiffness in the 
beam-column joint. 

Keywords: energy dissipation devices, structural fuses, reinforced concrete 
structures, pushover analysis, time-history analysis. 


1. INTRODUCCION 

Los ingenieros y academicos utilizan modelos matematicos para 
representar las caracterfsticas de interes de una estructura en cuestion y de esta 
manera establecen los parametros necesarios para su analisis y diseno. Con el 
avance tecnologico y computacional, dichos modelos matematicos tienden a ser 
mas complejos, debido a la incorporacion de nuevos algoritmos que simulan el 
comportamiento de distintos tipos de materiales o elementos estructurales. 

Se han realizado diversos estudios para el modelado de estructuras con 
disipadores histereticos de energia; sin embargo, estos se limitan a modelos 
pianos de uno o multiples grados de libertad. Por ejemplo, se han evaluado 
formulaciones para el modelado numerico de elementos ADAS, basados en la 
variacion de la seccion transversal del disipador con idealizaciones geometricas, y 
se han comparado con resultados experimentales obteniendo buena correlacion 
entre ellos (Whittaker et al. 1989, Tena-Colunga 1997). No obstante, existen otras 
propuestas sencillas, como las basadas en discretizar a los disipadores por medio 
de dovelas, donde cada una de ellas tiene longitud constante pero ancho variable 
(Aguiar etal. 2016). 

Otros estudios proponen metodologias de diseno basadas en desempeno, 
para sistemas equivalentes de un grado de libertad aplicables a sistemas de 
multiples grados de libertad. Dentro de estos se puede citar el trabajo publicado 
por Ruiz y Badillo (2001), donde evaluan edificios reforzados con disipadores de 
energia; o bien, la investigation Nevada a cabo por Chen y colaboradores para el 
diseno sismico de pilas de puentes estructuradas con marcos de acero y 
contravientos restringidos contra pandeo. En estos estudios se modela al sistema 
disipador con un comportamiento bilineal y se hace uso de analisis dinamicos paso 
a paso, ya sea para determinar las distorsiones de entrepiso y la ductilidad 
maxima de los disipadores en el sistema equivalente de un grado de libertad (Ruiz 
y Badillo 2001), o bien para validar la exactitud del procedimiento de diseno (Chen 
etal. 2007). 

Con el interes de evaluar la eficiencia de nuevas configuraciones o de 
nuevos materiales para los disipadores por histeresis, se han realizado 
investigaciones que constan tanto de pruebas experimentales como de estudios 
parametricos; estos ultimos, con el objetivo de comprobar la eficacia de dichas 
propuestas. A modo de ejemplo, existen estudios donde se proponen nuevas 
conexiones en forma de angulo para elementos de madera usados en 
edificaciones europeas de pequena a mediana altura. Estas conexiones se 
comportan bajo el mismo principio que los disipadores ADAS y analfticamente se 
representaron con elementos tipo liga no lineales definidos acorde al modelo 
multilineal de Sivaselvan y Reinhorn (Latour y Rizzano 2015). Tambien se puede 
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mencionar el estudio llevado a cabo para crear disipadores con forma de reloj de 
arena con el uso de cobre Cl 1000, que es ductil, de baja fluencia, y altamente 
resistente a la corrosion (De la Llera et al. 2004). A partir de ensayes 
experimentales de siete pares de disipadores, se identificaron los parametros que 
controlan su comportamiento ciclico y se propuso un modelo simplificado 
dependientes de su relacion de aspecto y de su ductilidad. Este modelo propuesto 
se uso para modelar a los disipadores de cobre (definidos por una unica 
deformacion transversal), en marcos pianos de concreto reforzado de 6, 12 y 25 
niveles con tres diferentes configuraciones estructurales (De la Llera et al. 2004). 

La calibration matematica entre resultados de marcos pianos con 
dispositivos ADAS, obtenidos a partir de programas de investigation y comerciales 
de analisis estructural no es un tema nuevo. Dentro de este ambito, se encuentra 
la comparativa hecha por Chaleo (2014) con elementos finitos elaborados en 
MATLAB y evaluados con dos ejemplos elaborados en el programa SAP2000. 
Chaleo establece una idealization con elementos finitos para el dispositivo ADAS 
con elementos de liga para representar su comportamiento histeretico bilineal. 
Cada propiedad del elemento ADAS idealizado se puede atribuir a un solo 
desplazamiento interno, es decir, su deformacion no afecta al comportamiento de 
otros desplazamientos internos. Por tal motivo, se aprecian diferencias en la 
comparacion de las respuestas estructurales mientras se utilizan un mayor numero 
de dispositivos ADAS y propone que se obtendrfa una mejor correlation de los 
resultados si se hubiese modelado a los disipadores con base en el modelo Bouc- 
Wen. 


En trabajos previos de los autores (Tena y Nangullasmu 2013; 
Nangullasmu y Tena 2015), se evaluo el comportamiento de marcos pianos con 
disipadores de energia a partir de una metodologia basada en balances de rigidez. 
Estos resultados se obtuvieron de analisis estaticos no lineales efectuados en el 
programa DRAIN 2DX (Prakash et al. 1992), modelando a los disipadores como 
elementos viga-columna. Los resultados de dicho programa son altamente 
confiables y se han calibrado en un sinnumero de investigaciones anteriores a 
esta, pero cuenta con varias limitantes, entre ellas la ausencia de interfaces 
graficas y el realizar los analisis en dos dimensiones. Ademas, para visualizar el 
desempeno de los marcos paralelos al del objeto en estudio, se tienen que hacer 
uso de otras tecnicas de modelado. Con el objetivo de realizar modelos mas 
complejos en 3D, que tomen en cuenta efectos ortogonales, se comparan los 
resultados de una serie de modelos realizados en DRAIN 2DX con el programa 
MIDAS Gen. Estos resultados se obtuvieron a partir de analisis estaticos no 
lineales y de dinamicos paso a paso. 


2. DESCRIPCION DE LOS MODELOS EVALUADOS 
2.1 Configuracion geometrica 

Se analizaron marcos pianos, marcos pianos acoplados y modelos 
tridimensionales de concreto reforzado con disipadores de energia soportados 
mediante un sistema de contraventeo de acero tipo chevron. La configuracion 
geometrica de los modelos en estudio corresponde a edificios regulares; se 
analizaron modelos de cinco y 20 niveles de marcos pianos formados por cuatro 
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crujfas de 800 cm de claro, alturas de entrepiso de 400 cm y angulo de inclination 
de contravientos 0=45° (Figura 1). Los modelos de marcos pianos acoplados 
tienen la misma geometrfa que los marcos pianos anteriormente descritos; se 
modelaron cuatro marcos en la direccion “Y’ con cuatro crujfas (Figura 2) y cinco 
en la direccion “/ con tres crujfas (Figura 3), acoplandolos mediante vigas enlace 
axial e infinitamente rfgidas para cada direccion en estudio, y de esta manera 
visualizar el efecto de los marcos paralelos en cuestion para compararlo con los 
resultados obtenidos del modelo tridimensional de cinco niveles (Figura 4). 



a) Modelo de cinco niveles 

Figura 1 Modelo de marcos pianos 
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Figura 2 Modelo de marcos acoplados direccion V 



Figura 3 Modelo de marcos acoplados direccion “/ (dimensiones en cm) 
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Figura 4 Modelo tridimensional 

Los cambios de colores en las Figuras 1 a 4 representan 
esquematicamente los cambios de secciones transversales conforme aumenta la 
altura de los modelos, es decir, los colores comunes en uno y otro modelo no 
corresponden a una misma seccion (por ejemplo, las vigas amarillas del modelo 
de cinco niveles son distintas a las vigas amarillas del modelo de 20 niveles). En 
ninguno de los modelos existe cambio de secciones de vigas, columnas y 
contravientos en el mismo nivel de entrepiso, con el objetivo de mitigar el efecto de 
piso suave por el cambio abrupto de rigidez y resistencia. 

2.2 Configu radon y modelado estructural 

Para la comparacion de resultados, usando dos programas distintos de 
analisis estructural, se modelaron a los elementos estructurales siguiendo la 
misma hipotesis de comportamiento para cada sistema estructural. Los elementos 
de contraventeo se modelaron para que trabajasen exclusivamente en direccion 
axial. Vigas y columnas se modelaron tomando en cuenta efectos de rigidez para 
tension/compresion, cortantes y momentos flexionantes. 

El modelado de los disipadores de energia se puede hacer de diversas 
maneras dependiendo del programa de analisis no lineal utilizado y de los 
elementos disponibles en sus librerias (Tena 2000). Conforme los programas de 
analisis estructural han ido evolucionando, dichos elementos disponibles se han 
ampliado con el paso del tiempo. Por ejemplo, era comun a principios de los anos 
noventa modelar los disipadores de manera simplificada con elementos axiales 
equivalentes, modelado que para un montado en contraventeo chevron es 
razonable, ya que los elementos equivalentes trabajan en la misma direccion que 
lo hacen los dispositivos, aunque se obtienen imprecisiones debido a la interaccion 
artificial que se establece entre los nudos superiores de los marcos de cada 
entrepiso con los disipadores, lo que en realidad no sucede al estar los 
disipadores conectados exclusivamente en la parte intermedia de las vigas (Tena 
2000). 


Despues, al desarrollarse elementos tipo liga (“link”) de cero longitud que 
se discretizan en un nodo, estos se han utilizado para modelar disipadores de 
energia (por ejemplo, Christopoulos y Filiatrault 2006). La ventaja es su sencillez 
de aplicacion. La desventaja en un montado con contravientos chevron es que no 
modela la longitud real de los disipadores, generalmente se consideran que solo 
trabajan a cortante (se desprecia la flexion y la carga axial, o si se consideran, 
estan totalmente desacopladas por ser modelados por medio de resortes), y no se 
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considers la rotacion que puede experimentar estos en su apoyo inferior con el 
sistema de contravientos, que modifica el angulo real de inclination de estos (que 
ya se ha demostrado impactan a la respuesta real del marco) y, por ello, se tienen 
imprecisiones que no han sido cuantificadas ni reportadas aun en los elementos 
mecanicos de los contravientos y de las vigas que se conectan directamente con 
los disipadores. 

Siempre ha existido la posibilidad de modelar con mayor rigor a los 
disipadores (Tena 2000), utilizando elementos tipo viga-columna (“beam-column”) 
que permiten considerar su longitud, ubicacion y condiciones de continuidad 
“reales”, asi como considerar que pueden trabajar a cortante, flexion y axialmente 
al mismo tiempo. El esfuerzo adicional para hacerlo no es muy grande, es incluir 
los nodos extras requeridos para modelar con precision la geometrfa del disipador, 
y los resultados que se obtienen correlacionan muy bien con los que se han 
reportado en experimentos que se han hecho en mesa vibradora (Tena-Colunga 
1997). Sin embargo, es un hecho que la popularidad de los elementos tipo liga 
para modelar disipadores ha aumentado en anos recientes. 

Asf, en este estudio los disipadores de energfa se modelaron en el 
programa DRAIN 2DX como elementos tipo viga-columna, mientras que en MIDAS 
Gen (Midas 2014) se modelaron de dos maneras distintas. La primera de ellas, 
como elementos tipo viga-columna y su formulacion esta basada en la teoria de 
vigas de Timoshenko (Figura 5a). La segunda, como elementos tipo liga, los 
cuales son resortes que individualmente representan una deformacion axial, dos 
deformaciones por cortante, una deformacion por torsion y dos deformaciones por 
flexion (Figura 5b). Dependiendo del metodo de aplicacion del analisis, puede ser 
del tipo elemento o del tipo fuerza. Los elementos de liga del tipo elemento reflejan 
directamente el comportamiento no lineal del elemento restableciendo la matriz de 
rigidez (analisis estaticos no lineales). Cuando se utilizan los elementos de liga del 
tipo fuerza (analisis dinamicos paso a paso), estos reflejan la no linealidad de 
manera indirecta convirtiendo las fuerzas de los elementos, calculadas a partir de 
las propiedades no lineales, a fuerzas externas (Midas 2014). 




Figura 5 Modelado de elementos estructurales en MIDAS Gen (adaptado de Midas 2014) 

Otra diferencia importante del modelado entre los programas utilizados 
radica en el punto de insercion geometrico de las vigas. Mientras que en DRAIN- 
2DX este punto de insercion viene definido de manera predeterminada y se puede 
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ajustar con la ubicacion de las zonas rfgidas de los extremos de las columnas, en 
MIDAS Gen se puede elegir directamente entre insertar la viga en el centroide 
geometrico de la seccion o en el punto superior medio para alinearlo con la altura 
del entrepiso. Por este motivo y como primera instancia, se evaluaron cinco 
versiones del modelo de marco piano de cinco niveles a partir de los resultados de 
analisis estaticos no lineales. Cuatro de estas versiones se modelaron en MIDAS 
Gen en funcion de dos parametros; el primero, la forma de modelar el disipador 
(elementos viga-columna o elementos de tipo liga); la segunda, la variation del 
punto de insercion de las vigas. Para todos los modelos las condiciones de apoyo 
para las columnas del primer nivel son empotradas y se agrego una condition de 
frontera a giro en el extremo inferior de todos los disipadores. 

La finalidad de la evaluation de los marcos pianos de cinco niveles fue 
valorar la divergencia de resultados relacionados al modelado entre ambos 
programas. En los marcos pianos de 20 niveles se aprecian las diferencias en el 
modelado del disipador para los dos tipos de analisis efectuados en este trabajo 
(estatico no lineal y dinamico paso a paso). Finalmente, y con base en la 
comparativa de los resultados de marcos pianos, se evalua el impacto de la rigidez 
del nudo viga-columna en la rigidez global a partir de las diferencias observadas 
entre los marcos pianos acoplados y el modelo tridimensional llevados a cabo en 
DRAIN-2DX y MIDAS Gen respectivamente. 

Las siglas utilizadas para los modelos hacen referenda al numero de 
niveles usados (M5), al tipo de elemento usado para los disipadores (B para 
elementos viga-columna y L para elementos de tipo liga), y al punto de insercion 
de las vigas en el modelo (C para el centroide y T para la parte superior de la 
seccion transversal). A manera de ejemplo, el modelo M5BC representa al marco 
de cinco niveles, con disipadores modelados con elementos viga-columna y las 
vigas colocadas en el centroide de la seccion transversal. 

Los modelos evaluados en el presente trabajo se disenaron con base en 
relaciones de rigidez entre el sistema del marco y el sistema contraviento-disipador 
(a) y entre el contraviento y el disipador (p). Ademas de cuidar que se cumplieran 
los balances de rigidez antes mencionados, se revisaron que todos los elementos 
cumplieran con las demandas de resistencia impuestas en el proceso de 
prediseno (Tena y Nangullasmu 2013). Los elementos de concreto reforzado se 
disenaron acorde a los requisitos de marcos no ductiles de las Normas Tecnicas 
Complementarias (NTCC-04), anadiendo al proceso de diseno los requisitos de 
confinamiento ductil referidas a las secciones 7.2.3 y 7.3.4, relacionadas 
principalmente a la separation de estribos en vigas y columnas respectivamente. 

Se presenta en la Tabla 1 las secciones utilizadas para cada uno de los 
sistemas involucrados para los modelos de marcos pianos acoplados de cinco 
niveles y para el modelo tridimensional. Los balances de rigidez indican que el 
sistema del marco es mas flexible que el sistema estructural global (a=0.25), y que 
el disipador tiene la mitad de rigidez que el contraviento (P=0.50). De manera 
esquematica, en la Figura 6 se aprecian los armados longitudinales y 
transversales de las secciones de concreto. El diametro de los estribos utilizados 
se distinguen con colores diferentes; rojo, para estribos No. 2.5 (diametro de 5/16” 
o 7.9 mm); azul, para estribos No. 3 (diametro de 3/8” o 9.5 mm); verde, para 
estribos No. 4 (diametro de 1/2” o 12.7 mm). 
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Se aprecian de la tabla y figura antes mencionadas, que para vigas y 
columnas las secciones transversales variaron conforme a la altura de los 
modelos. Ademas, se decidio hacer diferencias en el armado longitudinal para las 
vigas acorde a la direccion de analisis, y en las columnas se reforzaron en mayor 
medida aquellas donde se descargaba la componente axial de los contravientos. 
Como se aprecia de las Figuras 2 y 3, la rigidez de los marcos de concreto es 
diferente para cada direccion de analisis, debido al numero de crujias que 
presentan cada uno. Por tal motivo, se decidio controlar los balances de rigidez a 
partir del espesor de los contravientos de acero utilizados y de la rigidez elastica 
de los disipadores, para poder satisfacer de manera practica los balances “a” y “/3” 
propuestos. 


Tabla 1 Modelo de 5 niveles, a=0.25, p=0.50 (unidades: t-cm) 



Entrepiso 

b 

h 

Lecho Superior 

Lecho Inferior 

Est. Extremo 

Est. Centro 

V) 

Vlx 1-3 

35 

55 

6 Vs No. 6 

4 Vs No. 6 

4 No. 2.5 @ 10 

4 No. 2.5(g) 15 

(0 

CT) 

V2x 4-5 

30 

50 

5 Vs No. 6 

4 Vs No. 6 

3 No. 2.5 @ 10 

4 No. 2.5(g) 15 

> 

Vly 1-3 

35 

55 

7Vs No.6 + 4Vs No.5 

6 Vs No.6 + 2Vs No.5 

4 No. 2.5 @ 10 

4 No. 2.5(g) 15 


V2y 4-5 

30 

50 

6Vs No.6 + 5Vs No.5 

6 Vs No.6 + 2Vs No.5 

4 No. 2.5 @ 10 

4 No. 2.5 @ 15 


Entrepiso 

b 

h 

Refuerzo Long. 

Cuantia 

Est. Extremo 

Est. Centro 

(0 

c 

Cl Ext 1-3 

60 

60 

48 Vs No. 6 

3.80% 

7 No. 3 @ 10 

7 No. 3(g) 15 

E 

3 

C2 Ext 4-5 

55 

55 

40 Vs No. 6 

3.77% 

6 No. 3 @ 10 

6 No. 3(g) 15 

O 

o 

Cl Int 1-3 

60 

60 

20 Vs No. 6 

1.58% 

4 No. 3 @ 8 

4 No. 3(g) 15 


C2 Int 4-5 

55 

55 

16 Vs No. 6 

1.51% 

3 No. 4(g) 10 

3 No. 4(g) 15 

o 

Entrepiso 

b 

h 

t 




c 

Cl) 

CVIx 1-2 

30.48 

30.48 

1.91 




> 

CO 

CV2x 3-5 

25.4 

25.4 

1.91 




c 

o 

CVIy 1-2 

30.48 

30.48 

1.59 




o 

CV2y 3-5 

25.4 

25.4 

1.27 





Entrepiso 

Fy 

k 2 

KeLASTICA 




3— 

o 

■a 

Dlx 1-2 

65.41 

0.05 

251.34 




CO 

Q. 

D2x 3-5 

51.77 

0.05 

206.66 




'a> 

a 

Dly 1-2 

65.41 

0.05 

211.78 





D2y 3-5 

51.77 

0.05 

141.50 





V1X V2X V1Y 


V2Y 



Figura 6 Armado de vigas y columnas para el modelo de 5 niveles con a=0.25, p=0.50 
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2.3 Diferencias en la introduction de datos 

El programa MIDAS Gen cuenta con dos modalidades para la introduction 
de los diagramas momento-rotacion y de los diagramas de interaccion, ambos 
indispensables para los analisis estaticos no lineales. El primero de ellos, es la 
forma automatica a partir del armado propuesto por el usuario en las secciones de 
concreto utilizadas. La segunda es introduciendo de manera manual los datos 
para crear los diagramas antes mencionados, uno por cada tipo de section a 
considerar. En la Figura 7 se presentan los diagramas momento-rotacion 
obtenidos a partir del calculo automatico en MIDAS Gen (Ifnea continua) y los 
calculados de manera externa. Se aprecia de esta, que MIDAS Gen no cuenta aun 
con una herramienta necesaria para el calculo de los diagramas que tome en 
cuenta el confinamiento del concreto, causado por la separation de los estribos. 



Figura 7 Diagrama momento-rotacion para vigas 

Para cada una de los modelos antes mencionados, se introdujeron de 
manera manual los datos de los diagramas momento-curvatura y de los diagramas 
de interaccion calculados previamente con el programa BIAX (Wallace y Moehle 
1989), con el objetivo de considerar el confinamiento del concreto reforzado y 
obtener una mejor correlation entre los resultados de ambos programas. La 
manera de introducir los parametros requeridos para el analisis es sustancialmente 
diferente para cada uno de los programas utilizados; mientras DRAIN-2DX 
requiere de un formato de texto, MIDAS Gen cuenta con una interfaz grafica la 
cual permite una visualization previa de los datos introducidos. 

Para los diagramas de momento-rotacion se realizo una idealizacion 
bilineal (Figura 8), asegurando que la cantidad de energia entre el diagrama 
obtenido por el programa BIAX fuera equivalente al del idealizado. Ambos 
programas necesitan el momento de fluencia del elemento y su pendiente post- 
fluencia. La diferencia principal radica en que DRAIN-2DX no delimita el momento 
ultimo que resiste el elemento, mientras que en MIDAS Gen se puede plantear 
dicho momento. De esta manera, MIDAS Gen determina el punto de quiebre de la 
curva de capacidad del sistema global y, con ello, el punto ultimo del analisis 
donde la estructura se vuelve inestable. 
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Viga 2 



Fig ura 8 Idealizacion del diagrama momento-rotacion para viga 

Para los diagramas de momento-rotacion; en DRAIN-2DX, se hacen 
idealizaciones entre los puntos de la falla balanceada y los de compresion axial 
pura y momento en ausencia de carga axial (linea roja). Por otro lado, en MIDAS 
Gen se permite la introduccion de un mayor numero de datos del diagrama de 
interaccion, obteniendo un suavizado del diagrama que se asemeja con mayor 
precision al del calculado con BIAX (Figura 9). 


Columna 3 



Figura 9 Diagrama de interaccion momento-carga axial para columna 


3. ANALISIS ESTATICOS NO LINEALES 

Con base en analisis estaticos no lineales bajo cargas laterales monotonas 
crecientes, efectuados con el programa DRAIN-2DX y MIDAS Gen, se presentan 
resultados caracteristicos para los diversos modelos evaluados. Se presentan 
comparativas entre mapeo de fluencias y curvas de cortante basal normalizado vs 
distorsion global de azotea. Ademas, se reportan en tablas las diferencias en los 
parametros sismicos de diseno obtenidas a partir de las curvas mencionadas con 
los diferentes parametros de modelado tornados en cuenta. 
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3.1 Marcos pianos de cinco niveles 

Con el objetivo de evaluar las distintas formas de modelado de los 
disipadores en MIDAS Gen y de observar la posible divergencia de resultados 
ocasionado por el punto de insercion de las vigas, se evaluan en primera instancia 
marcos pianos de cinco niveles con las caracterfsticas de la Figura la. Para los 
modelos en MIDAS Gen donde se utilizaron elementos liga, estos fueron del tipo 
elemento (M5LC y M5LT). En la Figura 10 se presentan las curvas cortante basal 
normalizado-distorsion global de azotea; con Ifnea solida roja, se indican los 
resultados obtenidos con el programa DRAIN-2DX; con Ifneas formadas por 
guiones, los obtenidos para marcos con disipadores modelados con elementos 
viga-columna (M5BC y M5BT); y con Ifneas formadas por guiones y puntos, los 
calculados para marcos con disipadores modelados con elementos liga. 


Modelo de 5 niveles 



Distorsion (%) 


Figura 10 Curva cortante basal normalizado-distorsion global de azotea con a=0.25, p=0.50 

De los resultados mostrados en la Figura 10, y como era de esperarse, se 
aprecia que el modelar a los disipadores con elementos viga-columna en MIDAS 
Gen hace que la rigidez elastica del modelo global coincida exactamente con los 
resultados obtenidos por el DRAIN 2DX puesto que se empleo el mismo principio 
de modelado para los disipadores entre ambos programas. Sin embargo, se 
observa una variation perceptible en el intervalo no lineal, sobre todo, si se evalua 
el cambiar el punto de insercion de las vigas en MIDAS Gen, obteniendo mayor 
magnitud de sobrerresistencia cuando se modela a la viga para que coincida con 
la altura de entrepiso (M5BT y M5LT), sin importar el tipo de elemento que se 
utilice para modelar a los disipadores. Por otro lado, modelar a los disipadores con 
elementos liga en MIDAS Gen, hace ligeramente mas rfgida a la estructura global, 
debido a que no se modela la rotation que experimenta el disipador con respecto 
al contraviento de soporte. 

Las diferencias en los resultados debido al cambio en el punto de insercion 
de las vigas, tienen su fundamento en la forma de modelar las zonas rfgidas en 
MIDAS Gen. La longitud de rigidez infinita en las columnas varfa segun el punto de 
insercion, es decir, existen diferencias si hace coincidir el nivel de entrepiso con la 
parte superior de la viga o con el centroide geometrico de esta (Figura 11a). 
Cuando se toma en cuenta la primera suposicion, la magnitud de la longitud de 
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rigidez infinita se calcula en el extremo superior de la columna (Figura 11b). Si se 
toma en cuenta la segunda hipotesis, la longitud de rigidez infinita se calcula para 
ambos extremos de la misma. Ademas, al variar dicho punto de insertion, se 
modifica en automatico la rigidez del elemento y la distribution de carga, es decir, 
se crea una pequena excentricidad, lo que genera fuerzas cortantes y momentos 
flexionantes adicionales. Esta variacion no afecta en el calculo de deformaciones 
axiales o torsionales (Midas 2014). 


/ Zona Rigida 


I 


Longitud de rigidez infinita 


Zona Rigida 


a) offset centro-centro 


•a 



Longitud de rigidez infinita 



b) offset: centro-superior 


Figura 11 Modelado de zonas rigidas en Midas Gen (adaptado de Midas 2014) 


Evaluando los parametros sfsmicos de ductilidad (Q) y sobrerresistencia 
(R) reportados en la Tabla 2, los cuales se obtuvieron a partir de cada una de las 
curvas mostradas en la Figura 10, se comparan los diferentes resultados para 
cada programa en cuestion. Se aprecia que el modelar los disipadores bajo las 
mismas caracterfsticas entre los dos programas (elementos viga-columna), las 
diferencias no son mayores al 7% (sombreados en gris). Inclusive, para los 
modelos donde se presentaba una diferencia en la rigidez elastica por modelar a 
los disipadores de energia con elementos liga, las discrepancias de resultados no 
fueron mayores al 8%. El modelo que presento mayor variacion con respecto a los 
otros fue el M5BT modelado en MIDAS Gen. El uso de elementos viga-columna, 
aunado con alinear la viga con la altura de entrepiso, permitio al modelo 
desarrollar un mayor nivel de sobrerresistencia de mas del 14%, comparandolo 
con el modelo elaborado en el DRAIN-2DX. 


Tabla 2 Comparativa de resultados entre DRAIN-2DX y MIDAS Gen, marcos pianos cinco 
niveles 


PROGRAMA 

MODELOS 

RESULTADOS 

Q R T (s) 

% PORCENTAJES 

Q R T(s) 

DRAIN-2DX 

M5BC 

5.89 

1.73 

0.634 





M5BT 

5.62 

1.97 

0.632 

95.3% 

114.1% 

99.7% 

MIDAS 

M5BC 

5.62 

1.83 

0.633 

95.3% 

106.2% 

99.8% 

M5LT 

5.48 

1.80 

0.591 

92.9% 

103.9% 

93.2% 


M5LC 

5.55 

1.61 

0.598 

94.1% 

92.9% 

94.3% 
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3.2 Marcos pianos de 20 niveles 

Con base en los resultados evaluados para los modelos descritos de los 
marcos pianos de cinco niveles, se calibro el modelo de 20 niveles con balances 
de rigidez a=0.25 y 3=0.50 representative de la Figura 1b. Se decidio utilizar dicho 
modelo debido a que cuenta con un periodo cercano a 2 segundos (T=2.0 s), con 
lo que estarfa cercano al periodo dominante de la zona para el cual fue disenado 
(zona lacustre de la Ciudad de Mexico), y de esta manera, aprovechar el mismo 
para los analisis dinamicos paso a paso y hacer la comparativa de los ciclos de 
histeresis de los disipadores obtenidos con DRAIN-2DX y MIDAS Gen, que se 
presentan y discuten mas adelante. 

En la Figura 12 se grafican las curvas de cortante basal normalizado vs distorsion 
global de azotea de los resultados obtenidos con el programa DRAIN-2DX (Ifnea 
continua) y los modelados con MIDAS Gen (Ifnea discontinua). Se aprecia de esta, 
que existen diferencias en la rigidez elastica de las curvas, directamente 
relacionada por la forma de modelar a los disipadores de energfa. Para la figura 
antes mencionada, los modelos realizados en MIDAS Gen se usaron elementos 
liga para simular la participacion de los disipadores dado que, como se aprecia de 
la Tabla 2, las discrepancias de resultados no afectarfan de manera practica en la 
obtencion de parametros sfsmicos de diseno. Como se observa de la Tabla 3, las 
mayores diferencias se presentan en la obtencion del parametro Q cuando se usa 
como punto de insercion de la viga el centroide geometrico de la misma. Por otra 
parte, cuando se alinea la viga con la altura de entrepiso la mayor diferencia se 
presenta al obtener el parametro R. 


Modelo de 20 niveles 



Figura 12 Curva cortante basal normalizado-distorsion global de azotea con p=0.50, a=0.25 

Tabla 3 Comparativa de resultados entre DRAIN-2DX y MIDAS Gen, marco piano de 20 
niveles 


PROGRAMA 


DRAIN-2DX 

MIDAS 


MODELOS 


RESULTADOS OBTENIDOS 

Q R T (s) 


% PORCENTAJES 

Q R T 


M20BC 

3.95 

1.66 

1.992 




M20LC 

3.29 

1.61 

1.866 

83.3% 

96.6% 

93.7% 

M20LT 

3.62 

1.83 

1.864 

91.6% 

110.3% 

93.6% 
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En la Figura 13 se presenta el mapeo de fluencias inelasticas para el 
marco piano de 20 niveles obtenidas con MIDAS Gen y DRAIN 2DX. Comparando 
la magnitud de las ductilidades desarrolladas por los disipadores con MIDAS Gen 
(Figura 13a), contra las obtenidas de los resultados del DRAIN-2DX (Figura 13c), 
son bastante similares refiriendonos al codigo de colores, con exception de los 
primeros dos niveles del marco evaluado. Para el codigo de colores de las 
fluencias inelasticas en vigas y columnas, se aprecia que en ambos modelos 
tienen participacion no lineal incipiente las columnas de la periferia izquierda. Por 
otra parte, apenas es perceptible el cambio de colores referentes a una 
participacion no lineal mayor en vigas de los entrepisos intermedios, debido al 
codigo de colores utilizados en MIDAS Gen, donde el tono de azul mas intenso 
representa una participacion lineal elastica del elemento. 



a) Disipadores en MIDAS b) Vigas y columnas en c) Post-proceso de DRAIN 
Gen MIDAS Gen 2DX 


Figura 13 Mapeo de fluencias inelasticas para marco piano de 20 niveles con p=0.50, 
0=0.25 


3.3 Marcos acoplados y modelos tridimensionales de cinco niveles 

Evaluadas las diferencias de modelado en marcos pianos entre MIDAS 
Gen y DRAIN-2DX, se realizaron modelos pianos acoplados de cinco niveles para 
comparar sus resultados con el modelo tridimensional. Para el modelo en MIDAS 
Gen se utilizaron elementos liga para simular a los disipadores y se eligio la parte 
media superior como punto de insertion de las vigas. En DRAIN-2DX se modelo a 
la estructura en dos dimensiones, por lo que se tuvo que realizar dos modelos por 
separado para ver el efecto de los marcos acoplados en la direccion “x” (Figura 2), 
y en la direccion “y” (Figura 3). La ventaja que presenta MIDAS Gen es que se 
modela a la estructura en tres dimensiones en un unico modelo, con lo que 
permite el analisis y la visualization de efectos ortogonales. 
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En la Figura 14 se presentan las curvas cortante basal normalizado- 
distorsion global de azotea para los modelos de cinco niveles antes mencionados. 
Se aprecia de esta dos figuras, que corresponden a la direccion “x” y “y”, 
respectivamente. En las curvas analizadas en la direccion “x”, se aprecia ademas 
de los resultados del modelo tridimensional en MIDAS Gen y del modelo acoplado 
efectuado en DRAIN 2DX, una tercera curva con llnea discontinua que representa 
la curva de capacidad con las mismas caracteristicas geometricas y de resistencia 
pero con la configuration estructural representativa de la Figura la. 


Direccion X 


Direccion Y 




Distorsion (%) 


Distorsion (%) 


Figura 14 Curva cortante basal normalizado-distorsion global de azotea, modelo de cinco 
niveles 


De los resultados expuestos en la Figura 14 para ambas direcciones, se 
observa que existen pequenas diferencias en la rigidez de cada una de las curvas 
graficadas. Como se expuso con anterioridad, esto se debe al efecto del modelado 
de los disipadores y, para estos modelos en particular, principalmente a la rigidez 
de los nudos viga-columna en una modelacion tridimensional rigurosa como la 
hecha en Midas (modelo 3Dx5M y 3Dy5M) con respecto a la hecha en marcos 
pianos en el DRAIN-2DX. Ademas, se aprecia que la distorsion ultima que 
desarrolla el conjunto de marcos acoplados analizados en DRAIN-2DX es 
ligeramente mayor que los resultados evaluados en MIDAS Gen. 

En las Tablas 4 y 5 se resumen las diferencias en la obtencion de 
coeficientes sfsmicos de disefio “O' y “fl” a partir de las curvas de capacidad de la 
Figura 14, y de los periodos de vibrar de los modelos obtenidos en cada uno de 
los programas. Se aprecia que la diferencia entre los valores es, para efectos 
practicos, despreciables en ambas direcciones. Es decir, existe buena correlation 
entre el modelado en tres dimensiones y el modelado acoplado en dos 
direcciones. La mayor diferencia entre los resultados obtenidos se encuentra en la 
comparativa del coeficiente sismico de ductilidad “O’. Se observa de la Tabla 4 
que el modelar un solo marco con disipadores trae como resultado que la 
ductilidad global que desarrollarfa la estructura sea mucho mayor en comparacion 
que si se toman en cuenta en el modelado los demas marcos involucrados en la 
misma direccion, independientemente si se modela con una tecnica de 
acoplamiento o en un modelo riguroso tridimensional. 
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Tabla 4 Comparativa de resultados de modelos en 2D y 3D, direccion X 


PROGRAMA 

MODELOS 

Q 

DIRECCION X 

R 

T (s) 

% DIFERENCIAS 

Q R T (s) 

MIDAS 

3Dx5M 

5.24 

2.36 

0.6361 



DRAIN-2DX 

2D5D 

3D5DAx 

5.71 

5.09 

2.23 

2.40 

0.6391 

0.6713 

9.1% 

2.9% 

5.5% 0.5% 

1 .8% 5.5% 


Tabla 5 Comparativa de resultados de modelos en 2D y 3D, direccion Y 


PROGRAMA 

MODELOS 

Q 

DIRECCION Y 

R T(s) 

% DIFERENCIAS 

Q R T (s) 

MIDAS 

3Dy5M 

5.00 

2.78 

0.6935 



DRAIN-2DX 

2D5D 

3D5DAy 

5.20 

2.84 

0.7303 

4.1% 

2.2% 5.3% 


Las Figuras 15 a 22 corresponden al mapeo de fluencias inelasticas de 
todos los elementos involucrados en el sistema estructural. Como se ha observado 
en los mapeos para los marcos pianos (Figura 13), los contravientos se mantienen 
elasticos durante todo el evento sismico, dado que serviran de soporte a los 
disipadores y, por ende, estos no pueden fallar en ningun momento. En la Figura 
15 se observa de manera esquematica, el mapeo de rotaciones para el modelo 
efectuado en MIDAS Gen. Cabe senalar nuevamente que en las graficas de 
MIDAS Gen los cfrculos en azul marino representan respuestas elasticas (Figuras 
15, 17 y 21). 



Figura 15 Mapeo de fluencias inelasticas de MIDAS Gen 3D para modelo de cinco niveles 

Comparando las ductilidades de los disipadores obtenidas con MIDAS 
Gen en la direccion “x” (Figura 16) con las obtenidas de los resultados del DRAIN- 
2DX (Figuras 18 y 19), se aprecia una diferencia apenas perceptible en el 
disipador del cuarto entrepiso. Ademas, se observa que el patron de formacion de 
fluencias inelasticas en vigas y columnas son muy similares, y el codigo de color 
para todos los modelos es de un dano incipiente en vigas de los entrepisos uno al 
tres, y rotaciones en la base de las columnas interiores del primer nivel por efecto 
de la condition de apoyo como empotramiento perfecto. 

De manera similar que para los mapeos en la direccion “x”, para la 
direccion ortogonal a esta, se aprecian leves diferencias en el codigo de colores 
asociado a la ductilidad desarrollada por los disipadores del entrepiso tres (Figuras 
20 y 22). Por otra parte, la mayorfa de las columnas presentaron comportamiento 
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elastico, con excepcion de las columnas internas del primer entrepiso de los 
marcos sin disipadores (Figuras 21 y 22). Finalmente, las vigas de los dos ultimos 
entrepisos en ambas direcciones de analisis, no presentaron comportamiento 
inelastico. 



Figura 16 Mapeo de fluencia de disipadores en MIDAS Gen para modelo de cinco niveles, 
direccion X 



Figura 17 Mapeo de fluencias inelasticas MIDAS Gen para modelo de cinco niveles, 
direccion X 



Figura 18 Mapeo de fluencias inelasticas DRAIN-2DX para modelo de cinco niveles, 
direccion X 



Figura 19 Mapeo de fluencias inelasticas DRAIN-2DX para modelo de cinco niveles sin 
acoplar, direccion X 




Figura 20 Mapeo de fluencia de disipadores MIDAS Gen para modelo de cinco niveles, 
direccion Y 
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Figura 21 Mapeo de fluencias inelasticas MIDAS Gen para modelo de cinco niveles, 
direccion Y 
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Figura 22 Mapeo de fluencias inelasticas DRAIN-2DX para modelo de cinco niveles, 
direccion Y 


De manera general, se aprecia que la ductilidad del sistema global (Q) es 
similar para ambas direcciones evaluadas (0=5), para los modelos acoplados y 
tridimensionales. Por otra parte, existen mayores diferencias en los valores de 
sobrerresistencia obtenidos, y se aprecia que el valor de “ FT es menor en la 
direccion “x” que en la direccion debido a que en esta ultima se cuenta con 
mayor numero de crujias que influyen en una mejor distribution de la resistencia 
en los elementos involucrados. 


4. ANALISIS DINAMICOS PASO A PASO 
4.1 Filtrado de senales 

Con el objetivo de evaluar la participation de los disipadores en un analisis 
dinamico paso a paso con el programa DRAIN-2DX y compararlo con el programa 
MIDAS Gen, se procedio como primer paso a realizar un filtrado y correction de 
todos los registros de aceleracion utilizados en los modelos evaluados. El proceso 
de filtrado y correccion se realizo conforme lo descrito en Gomez et al. (2005), 
donde se establece lo siguiente: 

1. Correccion de linea base del registro en aceleracion, por posibles errores 
de medicion (correccion de linea base). 

2. Filtrado de la serial en aceleracion. Se uso un filtro pasa banda con una 
frecuencia corte pasa-alta de 0.05 Hz y una frecuencia de corte pasa-baja 
de 30 Hz. 

3. Integracion del registro en aceleraciones ya corregido por linea base y 
filtrado, para obtener la historia del registro en velocidad. 

4. Correccion de la linea base del registro en velocidad por errores en el 
proceso numerico de integracion. 

5. Integracion del registro corregido en velocidad (por linea base) para 
obtener el registro en desplazamiento 

En la Figura 23 se presenta de manera esquematica las diferencias de emplear la 
correccion del registro de aceleracion bajo las pautas antes descritas y lo que 
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sucederia si se hiciera caso omiso de ellas. Se aprecia perfectamente que el solo 
filtrar la serial no es suficiente, pero no aplicar la metodologla completa antes 
expuesta puede acarrear errores al momento de procesar los resultados de los 
analisis dinamicos paso a paso. Por ejemplo, puede ocasionar incrementos 
crecientes y poco reales en velocidades y aceleraciones, mientras que un filtrado 
excesivo eliminana informacion importante del sismo (Gomez etal. 2005). 


Registro de Desplazamiento 



Tiempo (s) 

Figura 23 Registro de desplazamiento obtenido a partir del registro de aceleracion de SCT- 
85, WE 

4.2 Modelo Berkeley 

Con la finalidad de realizar y evaluar modelos mas complejos que 
representen las condiciones geometricas, de resistencia y las propiedades 
dinamicas de modelos reales, se han calibrado los resultados de modelos 
experimentales con los obtenidos a partir de analisis dinamicos paso a paso 
creados en MIDAS Gen. Se comparan los resultados presentados por Whittaker et 
al. (1989), obtenidos a partir de pruebas experimentales de un modelo de tres 
niveles con disipadores de energi'a tipo ADAS, sujeto al registro de aceleracion “El 
Centro-33” con direccion norte-sur, ensayado en la universidad de Berkeley en 
1988 (Figura 24a). 



b) 

Figura 24 Modelo de Berkeley con disipadores tipo ADAS 



Se realizaron dos modelos anallticos comparativos en MIDAS Gen 
representatives de la Figura 24b. Estos se realizaron en tres dimensiones con la 
configuration geometrica y los pesos de entrepisos reportados por Whittaker et al. 
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(1989). Los elementos de vigas y columnas fueron modelados con elementos viga- 
columna, los contravientos con elementos que trabajan exclusivamente a carga 
axial y los disipadores con elementos liga del tipo fuerza. La diferencia entre los 
dos modelos evaluados se establece en las propiedades que se dieron de alta 
para los disipadores de energfa. Para el modelo denominado “MOD_EXP”, se 
utilizaron las propiedades de rigidez y de resistencia de los disipadores reportados 
en el modelo experimental; para el modelo denominado “MOD_Tena”, se uso la 
formulacion matematica para los disipadores tipo ADAS propuesta por Tena- 
Colunga (1997). Es importante mencionar que en este ultimo (“MOD_Tena”), se 
considero un modelo numerico acoplado del disipador, es decir, se considero la 
interaction de la rigidez de cada una de las direcciones tomadas en cuenta. Caso 
contrario del modelo “MOD_Exp”, donde se considero un modelo desacoplado del 
disipador, dado que unicamente se contaba con datos relativos de rigidez y de 
resistencia de la direction principal donde trabajaban. 

En la tabla 6 se resumen los periodos de vibrar de los modelos evaluados 
en esta section; se aprecia de esta que existe una mejor correlacion del primer 
modo de vibrar para el modelo analitico que tomo en cuenta directamente los 
datos experimentales, aunque para el segundo modo existen diferencias mayores 
para ambos modelos analfticos. En las Figuras 25 y 26 se presentan las graficas 
de desplazamientos maximos y de distorsiones respectivamente, suscitados por la 
excitation del registro de aceleracion “El Centro-33”. Se aprecia de la Figura 25 
que con la formulacion matematica propuesta por Tena-Colunga (1997), se 
obtiene una buena correlacion entre los desplazamientos obtenidos con el modelo 
matematico y los resultados experimentales. Sin embargo, se observa de la Figura 
26 que los resultados de los modelos analfticos, comparandolos con lo reportado 
por Whittaker et al. (1989), presentan una diferencia notoria en la forma del perfil 
de distorsiones; no obstante, para fines practicos los valores numericos son muy 
similares entre sf. 


Tabla 6 Comparativa de periodos de vibrar de Modelos Berkeley 


Modelo 

T 1 

t 2 

T 7 Error (%) 

T 2 Error (%) 

Experimental 

0.47 

0.17 



MOD EXP 

0.4694 

0.2343 

0.13% 

37.82% 

MOD_Tena 

0.4131 

0.2343 

12.11% 

37.82% 
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Figura 25 Desplazamientos maximos 


Figura 26 Distorsiones de entrepiso 
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En la Figura 27 se aprecian los ciclos de histeresis de los disipadores, 
reportados en la prueba experimental, para cada nivel de entrepiso. Se aprecia de 
esta que los dispositivos disipan mayor cantidad de energfa en el primer nivel, es 
decir, el potencial de dano del sismo se distribuye sobre la altura del edificio. 
Ademas, pueden desarrollar desplazamientos maximos de hasta 10 veces su 
desplazamiento de fluencia de manera estable y sin degradacion del elemento. En 
la Figura 28 se presentan los ciclos de histeresis desarrollados por los disipadores 
con la formulacion matematica de Tena-Colunga (1997). Se observa que para 
cada uno de los ciclos de histeresis mostrados, existe buena correlation de 
resultados con pequenas diferencias en la ductilidad ultima desarrollada por los 
disipadores, pero la forma del ciclo y la pendiente posterior a la fluencia son 
identicas. 
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Figura 27 Ciclos de histeresis de los disipadores de prueba experimental (Whittaker et al. 
1989) 
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Figura 28 Ciclos de histeresis de los disipadores modelados acorde a Tena-Colunga (1997) 

4.3 Marco piano de 20 niveles 

Se realizo un analisis dinamico paso a paso al marco piano de 20 niveles, 
representative de la Figura 1b, con balances de rigidez a=0.25 y 3=0.50. Se 
comparan los ciclos de histeresis de los disipadores obtenidos con los programas 
DRAIN-2DX y MIDAS Gen. El modelo de 20 niveles se sometio a una excitation 
sfsmica producida por el registro de SCT 1985 E-W, con lo que el periodo del 
modelo estarfa cercano al periodo dominante del suelo T s =2.0 s (zona lacustre de 
la Ciudad de Mexico), para el cual fue disenado dicho marco. 

En la Figura 29 se presentan los ciclos de histeresis de los disipadores 
para cada entrepiso del modelo en cuestion. En el eje “x” se tienen las distorsiones 
como porcentajes y en el eje “/ el cociente del cortante actuante entre el cortante 
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con el cual se disenaron los disipadores. Para cada una de las graficas, se aprecia 
en color rojo los resultados obtenidos con el programa DRAIN-2DX y en color azul 
los resultados obtenidos con el programa MIDAS Gen. 
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Figura 29 Ciclos de histeresis de los disipadores sujetos al registro de aceleracion SCT-85 
WE 


Como se menciono en parrafos anteriores, la forma de modelar a los 
disipadores histereticos de energfa entre ambos programas fue con distintos tipos 
de elementos. Esta primera diferencia en el modelado ocasiona un ligero cambio 
de rigidez de los disipadores, que se aprecia en la comparativa de los ciclos de 
histeresis de estos para cada entrepiso en estudio (Figura 29). Este efecto 
ocasiona a su vez que la rigidez elastica del sistema completo sea diferente, tal 
como se aprecia en la curva cortante basal normalizado-distorsion global de 
azotea de la Figura 12. Comparando los periodos obtenidos por los dos 
programas, se aprecia que el periodo del marco con disipadores modelado con 
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MIDAS Gen es menor que aquel que se obtiene con DRAIN-2DX (Tabla 3), por lo 
que se comprueba que los disipadores modelados con DRAIN-2DX tienden a ser 
mas flexibles, dado que incluyen la rotation del disipador y desarrollan mas ciclos 
de histeresis y una ductilidad ultima mayor, sobre todo a partir del disipador 
referente al cuarto entrepiso, tal como se aprecia en la Figura 29. 

4.4 Marcos acoplados y modelos tridimensionales de cinco niveles 

Para los analisis dinamicos paso a paso del modelo acoplado y 
tridimensional de cinco niveles, se evaluaron distintos registros de aceleracion 
sinteticos para la zona de transicion de la Ciudad de Mexico (Godinez 2005, Tena- 
Colunga et al. 2007), con el fin de hacer coincidir el periodo fundamental del 
modelo analitico con el periodo donde se presenta la maxima aceleracion de los 
registros y, de esta manera, que los disipadores desarrollaran su maxima 
capacidad. Con los registros de aceleracion corregidos y filtrados, se procedio a 
realizar los espectros de pseudoaceleracion mostrados en la Figura 30, para las 
direcciones Este-Oeste y Norte-Sur respectivamente. La linea en negro representa 
el periodo desarrollado por el modelo de cinco niveles con disipadores de energia 
para las direcciones correspondientes. 

Espectro de Pseudoaceleracion Este-Oeste Espectro de Pseudoaceleracion Norte-Sur 




Figura 30 Espectros de pseudoaceleracion para registros sinteticos de la zona de 
transicion de la Ciudad de Mexico generados por Godinez (2005) 

Se decidio utilizar los registros de aceleracion sinteticos A024 y GR27 
(Figura 31), dado que las ordenadas espectrales de ambas estan cercanas al 
periodo dominante de los modelos analiticos a evaluar. A partir de estos registros 
de aceleracion sinteticos se obtuvieron los ciclos de histeresis de los disipadores 
del modelo de marcos pianos acoplados de cinco niveles en cada direction de 
analisis, y de esta manera hacer la comparativa con el modelado riguroso 
tridimensional. 

En las Figuras 32 y 33 se presentan los ciclos de histeresis desarrollados por los 
modelos de cinco niveles con disipadores de energia al inducirle el registro de 
aceleracion sintetico A024. En color verde, se presentan los resultados obtenidos 
con el programa DRAIN-2DX para los marcos acoplados; en color rojo, los 
resultados de un solo marco con disipadores representative de la Figura la; y en 
color azul, los resultados que se obtuvieron con el programa MIDAS Gen. 
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Figura 31 Registros sinteticos de aceleracion generados por Godinez (2005) 
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Figura 32 Ciclos de histeresis para el acelerograma A024 de los disipadores de energfa, 
direccion X 

Comparando los resultados de los ciclos de histeresis desarrollados por 
los disipadores colocados en la direccion “x” (Figura 32), se aprecia que el hecho 
de considerar todos los marcos paralelos en la direccion de analisis permite que 
los disipadores desarrollen mas ciclos de histeresis y una ductilidad ultima mayor, 
que comparandolo con el modelo donde se considera un solo marco. Comparando 
los periodos dominantes en la direccion “x” de los modelos del marco piano y del 
modelo tridimensional (Tabla 4), se aprecia que estos son similares y, por ende, 
los ciclos de histeresis para ambos modelos son parecidos entre si en cada uno de 
los entrepisos evaluados (curvas en rojo y azul). 
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Figura 33 Ciclos de histeresis para el acelerograma A024 de los disipadores de energfa, 
direccion Y 


Evaluando los ciclos de histeresis desarrollados por los disipadores en la 
direccion “y", obtenidos por MIDAS Gen (curvas azules) y DRAIN-2DX (curvas 
verdes), se presenta la misma tendencia que los resultados expuestos para los 
mismos modelos en la direccion “x”. Se observa tambien que para los resultados 
obtenidos con el registro A024, los disipadores desarrollan mas ductilidad en la 
direccion que en la direccion ortogonal a esta. 


5. CONCLUSIONES 

Con el objetivo de evaluar el comportamiento de estructuras de concreto 
reforzado con disipadores histereticos de energfa que asemejen, en mayor 
medida, las caracteristicas de geometria y de resistencia de modelos reales, se 
realizaron modelos de marcos pianos, marcos pianos acoplados y modelos 
rigurosos tridimensionales con los programas DRAIN-2DX y MIDAS Gen. 

En el programa DRAIN-2DX se modelaron los disipadores exclusivamente 
como un elemento tipo viga-columna idealizando propiedades de rigidez y de 
resistencia, dado que este es el modelado predilecto del segundo autor. Sin 
embargo, cabe senalar que se pueden usar en este programa otros modelados, 
entre ellas modelar tambien a los disipadores con un elemento tipo liga (“link”), 
identificado en su libreria como “connection element”. En MIDAS Gen, que es un 
programa mas moderno que dispone de interfaces graficas, se modelaron a los 
disipadores como elementos tipo liga (“general link’) y tambien como elementos 
viga-columna, que a diferencia del elemento tipo liga, permite modelar la condition 
de apoyo que le dan los contravientos chevron, permitiendo una rotation 
parcialmente restringida. Cuando se modela de esta manera en Midas Gen, los 
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resultados obtenidos entre Midas Gen y DRAIN-2DX son bastante similares entre 
si, como es de esperarse, pues utilizan las mismas hipotesis. 

Evaluando y comparando las curvas de cortante basal normalizado- 
distorsion global de azotea y los mapeos de fluencia inelastica, obtenidos de los 
analisis estaticos no lineales, se observa una buena correlacion entre los 
resultados obtenidos con ambos programas. De manera practica, los coeficientes 
para diseno sfsmico relacionados a ductilidad y sobrerresistencia son similares con 
leves diferencias ocasionados, como se menciono en el parrafo anterior, por la 
manera de modelar a los disipadores entre los programas. El detalle mas relevante 
entre ambos programas es el hecho de que con MIDAS Gen, los analisis se 
detienen en automatico cuando los elementos han llegado a su resistencia ultima, 
dado que se pueden definir diagramas momento-curvatura donde se precisa el 
punto y valor de la curvatura de falla; mientras que con DRAIN-2DX es necesario 
la revision y el post-procesamiento de los datos obtenidos para evaluar el punto 
ultimo donde los resultados aun son realistas, dado que los diagramas momento 
curvatura que se calculan internamente en el programa a partir de las curvas de 
interaction propuestas son idealizados a partir de la fluencia y, por ello, no se 
definen la curvatura asociada a la falla. 

Comparando los ciclos de histeresis desarrollados por los disipadores de 
energia, a partir de diversos analisis dinamicos paso a paso, se observa que 
existen diferencias en la ductilidad ultima desarrollada por los disipadores. Esta 
diferencia se va acentuando conforme el modelo es mas esbelto, pero sobre todo, 
esta relacionado con la variation de los periodos obtenidos entre ambos 
programas. 

Finalmente, se observo que el modelado de la estructura tomando en 
cuenta todos los marcos actuantes en la direction de analisis (marcos acoplados), 
da como resultado una ductilidad global “O’ menor que si solo se considera un 
marco de interes. Para el caso de los analisis dinamicos paso a paso, se observan 
diferencias en la distorsion alcanzada en los ciclos de histeresis desarrollados por 
los modelos analizados en MIDAS Gen y en DRAIN-2DX, obteniendo mayor 
capacidad de deformation con este ultimo. Esto se debe principalmente a que en 
un modelado tridimensional riguroso, se presenta mayor rigidez de los nudos viga- 
columna (en Midas Gen), pues esa zona rfgida se modela tridimensionalmente, 
con lo cual impacta de manera significativa la rigidez global de la estructura con 
respecto a la que puede hacerse en marcos pianos (DRAIN-2DX), donde la zona 
rfgida se modela en el piano conforme a la teorfa general de elementos viga- 
columna (por ejemplo, Tena 2007). 
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